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Abstract. This paper is concerned with the asymptotic behavior of the free energy for a class 
of Hermitean random matrix models, with odd degree polynomial potential, in the large TV limit. It 
continues an investigation initiated and developed in a sequence of prior works whose ultimate aim is 
to reveal and understand, in a rigorous way, the deep connections between correlation functions for 
eigenvalues of these random matrix ensembles on the one hand and the enumerative interpretations 
of their matrix moments in terms of map combinatorics (a branch of graph theory) on the other. In 
doing this we make essential use of the link between the asymptotics of the random matrix partition 
function and orthogonal polynomials with exponential weight equal to the random matrix potential. 
Along the way we develop and analyze the continuum limits of both the Toda lattice equations and 
the difference string equations associated to these orthogonal polynomials. The former are found to 
have the structure of a hierarchy of near-conservation laws universal in the potential; the latter are a 
novel semi-classical extension of the traditional string equations. Our methods apply to regular maps 
of both even and odd valence, however we focus on the latter since that is the relevant case for this 
paper. These methods enable us to rigorously determine closed form expressions for the generating 
functions that enumerate trivalent maps, in general implicitly, but also explictly in a number of cases. 
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1. Introduction The general class of matrix ensembles we want to analyze 
has probability measures of the form 

dfx tj = — 7^ exp{-iV Tr[Vj(M,tj)]}dM, where (1.1) 

^(A; *i) = ^A 2 +i^ (1.2) 

defined on the space T-L n of nxn Hermitean matrices, M, and with N a positive 
parameter. The normalization factor Z^\tj), which serves to make \it a probability 
measure, is called the 'partition junction of this unitary ensemble. Previous works, 



[HI EH HH U\ j have focussed on the case of even j for which the measure (1.1) is indeed 
normalizable for tj >0. The case of odd j is more complicated; it is not clear prima 
facie how to initiate a rigorous analysis in this setting. 

Very recently, however, a generalization of the equilibrium measure (which governs 



the leading order behavior of the free energy associated to (1.1)) was developed and 
applied to this problem, pjj. It is based on a complex contour deformation of the 
variational problem for the leading order of the free energy that was motivated by new 
ideas in approximation theory related to complex Gaussian quadrature of integrals 
with high order stationary points [4]. 

This analysis shows that an equilibrium measure associated to the weight 
exp{— -/V (|A 2 +<2i/+iA 2,y+1 ) }, with dominant exponent odd, will exist. It is con- 
structed explicitly for the case of a cubic weight, u = l, in [3]. A detailed study of 
the explicit construction for general odd-dominant weights, as opposed to just the 
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existence argument which may be deduced from [4], will be taken up elsewhere. The 
boundaries of the support of this equilibrium measure are determined by the simul- 
taneous solutions of the two equations: 



HA) 



A y/{X-A)[\-B) 

\V'{\) 
a y/(X-A)(X-B) 



dX = 



dX = 2ni. 



(1.3) 
(1.4) 



One can compute these integrals which, in the cubic case, V(A) = |A + t$X , leads 
to a pair of equivalent algebraic equations determining A and B. 

I(A + B) + 3i 3 (^ 2 + ^45+^ 2 )=0 (1.5) 



■3t a 



16' 



-A 6 



1 A 2 B+-AB 2 + -B 3 
16 16 16 



(1.6) 



It is natural to make the following change of variables: zq = A (B — A) 2 and uq 



h(A + B). The corresponding algebraic equations for zq, uq are 

u a + 3t 3 (u 2 + 2z Q ) = Q 
ul + 2z + 3t 3 (u 3 ) + 6u z ) = 2. 



(1.7) 
(1.8) 



With the above notations, the interval of support of the equilibrium measure in 
the cubic case may be written as 

[A,B] = [uo - 2^, u + 2</zH] . 

The equilibrium measure has a variational characterization i4j and from the variational 
equations the measure can be explicitly determined to be 



^- (1 + 3i 3 (A + u )) X [A,B] (A) y/(X-A)(X-B). 



A minimal basis for the ideal of relations given by (1.7) and (1.8) is 

3t 3 ul +u + 6t 3 z = 
-6t 3 z u + (l-z ) = 0. 

It is straightforward to use (1.10 1 to eliminate zq in ( |1.11[ ) and get 

18*3^ + 9t 3 ul + M + 6t 3 = . 



The resultant of (1.7) and (1.8) eliminating uq is given by 



R(z ) -- 



3i 3 1 6*3^0 
-6t 3 z 1-Zq 
-6t 3 z 1-zq 

:3t 3 (72t 2 3 z%-z 2 + l)=0. 



(1.9) 



(1.10) 
(1.11) 



(1.12) 



(1.13) 



(1.14) 



We note that (1.14) has a form that is reminiscent of the implicit equation for zq that 
one has in the case of even weights jTUl [7] . 
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For general polynomial potentials V with even dominant power, it is possible to 
establish the following fundamental asymptotic expansion [5] , [TO] of the free energy 
associated to the partition function. More precisely, those papers consider potentials 
of the form 

V(\) = \\ 2 + Y.^ J i (1-15) 

with J = 2v. Introducing a renormalized partition function, which we refer to as a tau 
function, 

(1 ' 16) 

where t = (ti,.. .tj) £M. J , this expansion has the form 



logr^jv (t) = n 2 e (x, t) + e x (x,t) + ^e 2 (x, t) + - + -^^e g (x, t) + ... (1.17) 



as n,N— >oo with held fixed. Moreover, for T= (1 — e, 1 + e) X 

({\t\<6}n{tj>0}) for some e>0,S>0, 

(i) the expansion is uniformly valid on a compact subsets of T; 

(ii) e g {x,t) extends to be complex analytic in T c — {(x,t) £ C' 7+1 1 \x — 1| < e, |f| < 5}; 

(iii) the expansion may be differentiated term by term in (x,t) with uniform error 
estimates as in <J±J> 

The meaning of ^ is that for each g there is a constant, K g , depending only on T 
and g such that 



iogT^jv (i)-n 2 e (x,t) ^—^e g (x,t) 

n, y 



- n 2g 



for (x,t) in a compact subset of T. The estimates referred to in ( fiii] ) have a similar 
form with N and ej(x,t) replaced by their mixed derivatives (the same derivatives 
in each term) and with a possibly different constant. 

This result is based on the analysis of a Riemann-Hilbert problem (RHP) for 
orthogonal polynomials on K whose exponential weight is associated to the weight 
of the random matrix measure. This RHP was first introduced in |llj for studying 
the asymptotic behavior of random matrix partition functions. The relevant analysis 
of this RHP for the above result was carried out in [3] by the method of nonlinear 
steepest descent [5J. In particular, in it is shown that the constants K g are explicitly 
determinable in terms of Airy asymptotics stemming from the Airy parametrix that 
is used in the vicinity of the endpoints of the support of the equilibrium measure to 
explicitly solve the RHP. 

This result extends directly to the case of V with odd dominant power (i.e. with 
J = 2^ + 1) once one has the existence of the equilibrium measure (which is explicitly 



given by (1.9| in the cubic case). This involves studying the asymptotic behavior of 
the appropriate non-Hermitean orthogonal polynomials for the given odd weight (sec 
section |2j. More precisely, the Riemann-Hilbert analysis of [9] carries over mutatits 
mutandis to a Riemann-Hilbert problem for the non-Hermitean orthogonal polynomi- 
als with the principal difference being that the contour along which the jump matrices 
are originally defined is no longer the real axis but rather a deformed contour [H [3] . 
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Our principal interest in this paper is to better understand the analytical structure 



of the coefficients e g for potentials of the form ( 1.15 1 when J there is odd. For J even 
these coefficients provide a wealth of information about problems in combinatorial 
enumeration as well as about eigenvalue correlations for random matrices [7]- One 
expects to see similar connections in the case of odd J but this is much less developed. 



Despite the fact that we are focussed on potentials of the form (1.2) that only 
depend on a single tj, we will need to appeal to properties for other parameters 

as well, specifically t\ and x. That is because in characterizing the e g we will want 
to make use of differential relations in t\ and tj between these coefficients as well as 
certain rescalings of these variables in terms of x: 

s 1 =x~^tx (1.18) 
s^xi-Hj. (1.19) 



We will also be studying the tau-functions (1.16) as functions of lattice variables 
on the non-negative integers, indexed by n, which depend analytically on (x,t) as 
parameters. Certain logarithmic derivatives of the tau functions with respect to these 
parameters satisfy difference equations which, in this context, we refer to as differ- 
ence string equations. Furthermore they satisfy differential (in tj) - difference (in n) 
equations classically known as the Toda lattice equations. Of particular relevance for 
describing and analyzing the e g will be the continuum limits of the difference parts of 
all these equations. These involve, as independent variables, Sj and s± as well as a con- 
tinuous " spatial" variable w in terms of which the differencing in our string and Toda 
equations may be regarded as a discretization. At the final stage, after one has recur- 
sively solved the continuum limit hierarchies for e g (and various of their derivatives) as 
functions of (x,Si,Sj,w), the auxiliary variables are set to (x,si,w) = (1,0,1) to arrive 
at the desired closed formulae for e g (sj) =e g (tj). (Note that when x = l,Sj = tj.) 

The continuum limit of the difference string equations is a hierarchy of nonlinear 
odes in w while that of the Toda equations is a hierarchy of quasi-linear pdes in Sj 
and w where in both cases g indexes the respective hierarchy. One will have these 
hierarchies for each value of j. 

We take a moment here to briefly explain the connection of the expansion 



(1.17) to combinatorial enumeration that was alluded to earlier. The e g (tj) = 
e^Jx = l,t = (0,...,0,ij)J (we have set J~j here) are generating functions for the 
enumeration of j -regular maps. A map is an embedding of a labelled graph into a 
compact, oriented and connected surface X with the requirement that the comple- 
ment of the graph in X should be a disjoint union of simply connected open sets. 
More specifically the asymptotic expansion coefficient e g is a generating function for 
enumerating (topological) equivalence classes of maps on a Riemann surface of genus 
g ((/-maps) whose embedded graphs are j-regular: 

m>l 

in which each of the Taylor expansion coefficients K g ^ (m) is the number of g-maps 
with m j-valent vertices. Consequently the Taylor coefficients of e g {tj),Hi'p{m), are 
non-negative integers. 

The notion of g-maps was introduced by Tutte and his collaborators in the '60s 
[15) as a means to study the four color conjecture. However, this subject soon took 
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on a life of its own as a sub-topic of combinatorial graph theory. In the early '80s 
a group of physicists [5] discovered a profound connection between the enumerative 
problem for labelled g-maps and diagrammatic expansions of random matrix theory. 
That seminal work was the basis for bringing asymptotic analytical methods into the 
study of maps and other related combinatorial problems. 

The trivalent case of map enumeration (which corresponds to the random ma- 
trix ensemble with cubic weight) is of particular relevance for problems in discrete 
geometry since the corresponding maps are dual to triangulations which are stable 
discretizations of the associated Riemann surfaces. In particular, in Section |4j in or- 
der to explicitly solve the continuum difference string equations we will use the fact 
that the coefficient zq appearing in the equilibrium measure is the generating func- 
tion for enumerating labelled trivalent maps on a sphere connected to two marked 
univalent vertices. This is dual to the generating function that enumerates ordered 
triangulations of the sphere with two marked loops. One also expects zq to have a 
combinatorial interpretation analogous to the one it has in the cases of even valence 
which is as a generating function for the Catalan numbers (in the case of valence 4) 
and generalized Catalan numbers in the cases of higher even valence [10l [7J . That 
interesting topic will be taken up elsewhere. 

Recently, in [7J, closed form expressions for all of the e g (tj) were derived for the 
cases of even j. This was based on the continuum limit of Toda lattice equations, 



developed in 10J, that are closely related to the random matrix ensembles (1.1). In 
this paper we will derive the analogous continuum j-Toda equations for odd j as well 
as the related hierarchy of continuum difference string equations. From these we will 
derive explicit closed form expressions for the e g {t^,) for some low values of the genus 
g. For example we will find that 

l w , , 1 (z -l)(zg-6*o-3) 

60 (<3) = 2 log(zo) + 12 WTv, ' 

ei (*s) = -^lo g (§-f), (1.21) 

1 (z 2 -l) 3 (4z 4 -93z 2 -261) 
2{ 3) 960 (^-3) 5 

where zq is given by the boundary of the equilibrium measure discussed above and is 



implicitly related to by the polynomial equation (1.14) 



1 — 7 2 — 79t 2 ? 3 

We expect that our methods will ultimately enable one to derive closed form expres- 
sions for the e g (tj) for all odd j and all genus g. In [3] the Taylor coefficients of eo 
and ei are calculated for the cubic case (j — 3). These derivations were based on a 
different approach using the classical string equations. 

The outline of this paper is as follows. In section [2] we summarize the necessary 
background on non-Hermitean orthogonal polynomials that is the basis for the validity 
of the asymptotic expansions that we study as well as their continuum limits. This 
section also presents a path formulation for both the Toda lattice equations and the 
difference string equations associated to these orthogonal polynomials. The latter 
in particular represent a novel method for the study of random matrix continuum 
limits. The continuum limits themselves are derived in Section [3j for general odd 
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valence, at least for the leading order and higher order homogeneous terms. To help 
make this paper more self-contained, a preliminary subsection of Section [3] is included 
that summarizes prior results on which the work in this paper is based. This part 
also contains a new result: a description of the asymptotic structure of the diagonal 
recursion coefficients for the orthogonal polynomials. This result was not needed 
previously because these diagonal coefficients vanish in the case of even potentials. 
Although this result has a similar character to what had previously been found for the 
off-diagonal recursion coefficients, the derivation is technically more complicated due 



to the fact that the Hirota expression (3.f2) for the diagonal coefficients is given in 



terms of a leading order differential-difference operator rather than the pure second 
derivative (3.131 for the off-diagonal coefficients. In Section [4] we specialize to the 
trivalent case and derive the full Toda and difference string equations up to order 
g= 1. Moreover, we illustrate the use of these methods by explicitly solving the g = l 
difference string equations. Finally in Section [5] we derive a recursive method for 
expressing each generating function e g (t^,) in terms of zq and use this to establish the 



explicit formulae (1.21) 



2. The Role of Orthogonal Polynomials and their Asymptotics 

Let us recall the classical relation between orthogonal polynomials and the space of 
square-integrable functions on the real line, K, with respect to exponentially weighted 
measures. In particular, we want to focus attention on weights that correspond to the 



the Hilbert space H = L 



NVJXJ 



random matrix potentials, V(A), (1.2), that interest us here. To that end we consider 



of weighted square integrable functions. This 



space has a natural polynomial basis, {7r„(A)}, determined by the conditions that 



7r„(A) = A" + lower order terms 
n n (X)w m (X)e~ NVW dX = forn ^ m. 



For the construction of this basis and related details we refer the reader to [5] . 

With respect to this basis, the operator of multiplication by A is representable as 
a semi-infinite tri-diagonal matrix, 



£ = 



(a 1 
b\ ai 

V 



a 2 



\ 



J 



(2.1) 



C is commonly referred as the recursion operator for the orthogonal polynomi- 
als and its entries as recursion coefficients. We remark that often a basis of or- 
thonormal, rather than monic orthogonal, polynomials is used to make this repre- 
sentation. In that case the analogue of (2.1) is a symmetric tri-diagonal matrix. 



As long as the coefficients {b n } do not vanish, these two matrix representations 
can be related through conjugation by a semi-infinite diagonal matrix of the form 

diag(l,&r\(W~\(&iM 3 ) ,•••)■ 

Similarly, the operator of differentiation with respect to A, which is densely defined 
on H , has a semi-infinite matrix representation, T>, which we now determine. Observe 
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that 

' <(A)7r m (A)e- w(A) dA = for n < m; 
ir' n (\)ir m {\)e~ NVW d\ = N f Tr n (\)V'(\)Tr m (\)e- NVW d\ for n>m 



= N J 7 r n (A){A + ^'- 1 }^ m (A) e - A '^dA; (2.2) 

hence, 

V = N(C+jtO- 1 )_ (2.3) 

where the "minus" subscript denotes projection onto the strictly lower part of the 
matrix. 

From the canonical (Heisenberg) relation on H, one sees that 

[dx,X\ = l, 

where here A in the bracket and 1 on the right hand side are regarded as multiplication 
operators. Using this and orthogonality one has 



^ {{d x ,\]ir n (\)}ir m (\)e- NV ^d\ = K n 6 n 
with n n > 0; 

= J n n (\){l\~d x ]7r m (\)e~ NV ^}d\, 

where we note that under this transposition of the bracket within the inner product 
the order of composition of the operators has interchanged and the minus sign on the 
derivative comes from integrating by parts; 



= [n n (\)Y,[C,V]^ m n e (\)e- NV Wd\, 



by (2.2| and (2.3); 



It follows that 



■K n [C,V] 



[C,V] = L 



From this observation one deduces a fundamental relation among the recurrence 
coefficients, 



i 

N' 



-I. 



(2.4) 



The relations implicit in (2.4) have been referred to as string equations in the physics 



literature, but their origins go further back to the classical literature in approximation 
theory [T3]. In fact the relations that one has, row by row, in (2.4) are actually 
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successive differences of consecutive string equations in the usual sense. However, by 
continuing back to the first row one may recursively de-couple these differences to get 
the usual equations. To make this distinction clear we will refer to the row by row 



equations that one has directly from (2.4 1 as difference string equations. 



C depends smoothly on the coupling parameter tj in the potential V(X) (see |1.2[ ). 
The explicit dependence can be determined from the fact that multiplication by A 
commutes with differentiation by tj and the following consequence of the orthogonality 
relations: 



K(A))7r m (A)e- Ary(A) dA = iV / A J '7r„(A)7r m (A)e- w(A) dA, forn > m. 



d_ 
dt. 



This yields our second fundamental relation on the recurrence coefficients, 



1 

N 



(2.5) 



which is equivalent to the j equation of the semi-infinite Toda Lattice hierarchy. 

2.0.1. Odd Weights and Non-Hermitean Orthogonal Polynomials 

When j is odd, H as defined above ceases to be a finite measure space; however, 
by deforming the real axis to an appropriate complex contour T one can define a 
non-Hermitean analogue of orthogonal polynomials with respect to this contour and 
weight, [SIS]- These polynomials may not be defined for all values of n but asymptoti- 
cally they exist (i.e., for n > no — 1 for some sufficiently large integer no) @]. Thus one 
can work on the space H = L 2 (T,e^ NV ^) of weighted square integrable functions 
on the deformed contour T. Of course in doing this deformation one can no longer 
relate the construction of the orthogonal polynomials to an inner product on H as 
was done before (hence the nomenclature non-Hermitean). Instead one works with 
the non-degenerate complex-valued bilinear form that integration naturally gives us. 
One can then, as we will shortly see, define a basis of polynomials whose pairwise 
product integrates to zero if they are of different degree. One can still use this basis 
to represent the recurrence operators and related operators through the bilinear form. 
To that end let 

iJ={span ofl,A,...A" 0_1 ,7r„ (A),7r„ 0+ i(A),...} 
where 7r„ (A) is a monic polynomial of degree n such that 



= y 7r„(A)A fc e- w(A) dAforfc = 0,...,n-l. 



With respect to this basis, multiplication by A is represented as 



C = 



do a x ... a„ _i 





V 



o 



\ 



^ri + l a n + l 1 

t>n +2 O no+ 2 



(2.6) 



where 7T no _i(A) = a no _iA n ° 1 H \-ai\ + a . 
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One may apply standard methods of orthogonal polynomial theory |14| to the 
lower right semi-infinite block of this matrix 

/ Ann 1 \ 



L = 



"no 

Ka + l a no + l 1 

K +2 a n +2 



V 



(2.7) 



7 



In particular there is a unique semi-infinite lower unipotent matrix A such that 

C = A- 1 eA 

where 

/0 1 \ 
00 1 







V 



7 



(For a description of the construction of such a unipotent matrix we refer the reader 
to Proposition 1 of [8].) 

This is related to the Hankel matrix 

fco ci c 2 . . . \_ 

ci c 2 c 3 . . 
c 2 c 3 c 4 . . 



H = 



\: : : -.J 



where 



is the k moment of the measure, by 



c k = / X k e- NV ^dX 



ADA* = 



( c n n + l c n +2 Crio+3 
c n a +2 c n +3 c n +i 
c no+3 c n +4 c n +5 

V ; i ; 



with 



D = diag {d na+1 , d no+2 . . . } 
defH n 



where T-L n denotes the nxn principal sub-matrix of % whose determinant may be 
expressed as (see Szego's classical text [Hj), 



det%„ =n\Z { £ ) (ti,t 2v +i) 



Zfif {t2v+l)~- 





7r 6XP< 


—N 2 









N 



1 " 

y ^ Vj^m'M, *2i/+l)- 

-^^log|A m -A, 



m— 1 



cTA, (2.6 
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where V(A;ti, i 2 „+i) = \\ 2 + t 1 \ + t 2v+ i\ 2v+1 . 

Remark 2.1. ^4s mentioned in the introduction, we will sometimes need to extend 
the domain of the tau functions to include other parameters, such as t\, as we have 
done here. Doing this presents no difficulties in the prior constructions. 

The diagonal elements may in fact be expressed as 



where 



'n.N 



d n (0) 



'n-l,N 



J N 



J N 



J N 



(0,0) 



J N 



(0,0) 



(2.9) 
(2.10) 



which agrees with the definition of the tau function given in (1.16). The second 



equality follows by reducing the unitarily invariant matrix integrals in (2.10) to their 
diagonalizations which yields (2.9) [S]. This also provides the connection to eigenvalue 
correlations alluded to in the introduction. Tracing through these connections, from C 
to D, one may derive the basic identity relating the random matrix partition function 
to the recurrence coefficients, 



h 2 

°n,N ' 



T 2 T 2 



N 



&njv(0) 



(2.11) 



which is the basis for our analysis of continuum limits in the next section. With this, 



the fundamental relations (2.4) and (2.5) continue to hold in the non-Hermitean case 
for n sufficiently large. 



Remark 2.2. It needs to be noted that the bilinear form used to define orthogonal 
polynomials and recurrence coefficients in this section depends on the choice of the 
contour T and therefore so do these polynomials and coefficients. However, it does 
not affect the asymptotics of these objects. This is a consequence of the fact that 
outside the locus of support of the equilibrium measure, one has exponential decay of 
the asymptotics. The deformation of T away from K is taken in these exponentially 
decaying regimes. We refer the reader to I 111 where similar issues concerning 
non-Hermitean orthogonal polynomials and their asymptotics are discussed but for a 
different problem. 

Remark 2.3. The fact that the lower degree recurrence coefficients may not exist 
in the non-Hermitean case creates technical difficulties in deriving the usual string 
equations since, as was pointed out earlier, this derivation requires that one be able 
to recursively relate higher degree recurrence coefficients all the way back to degree 0. 
However, this issue poses no problems for the asymptotic difference string equations 
nor for the asymptotic Toda equations which are what will be used in this paper. 
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2.0. 2. Path Weights and Recurrence Coefficients In order to effectively 
utilize the relations (2.4 2.5) it will be essential to keep track of how the matrix 
entries of powers of the recurrence operator, £ J , depend on the original recurrence 
coefficients. That is best done via the combinatorics of weighted walks on the index 
lattice of the orthogonal polynomials. The relevant walks here are Motzkin paths 
which are walks, P, on Z which, at each step, can increase by 1, decrease by 1 or stay 
the same. Set 

V 3 (m\,rri2) = the set of all Motzkin paths of length j from mi to m^. (2-12) 

Then step weights, path weights and the (rai,m.2)-entry of C 3 are, respectively, given 
by 

{1 if the p th step moves from n to n + 1 on the lattice 
a n if the p th step stays at n 

b\ if the p th step moves from n to n — 1 

cj(p)= n ^ 

steps pep 

P£PJ(mi,m 2 ) 

2.1. Motzkin Representation of the Difference String equations 

The difference string equations are given (for the j-valent case) by (2.4): 

C,(£+jt&- 1 )_]=~I. (2.14) 

This leads to a pair of equations: 

• the (n + l,n) entry gives 

0^(a n+1 ~a n )(C+jtC^) n+ln + (C + jtC^) n+2n 

+ (2.15) 

• and the (n,n) entry gives 

Jf = ( £ +^ _1 )„ +1 ,„ - ^'-" C ' \n : • ( 2 - 16 ) 

Let us work this out in terms of Motzkin paths for the particular case of j = 3. The 
equations for the diagonal and subdiagonal equations reduce respectively to 



x 
n 



{C"n+l,n C"n,n-l) + 3t (£„ +1 „ C n _ 



— (a n+ i — a n ) + 3t£^ l+l n ) + (£ n +2,n — C n +i, n -i) + St (£ n+2> „ — C n+ljn _i) 

where we have used the relation x - 



N- 



Referring to (2.12), we see that the relevant path classes here are 

V 1 (n + l,rt) = a descent by one step 
V 1 (n + 2 } n) =thc empty set 

P 2 (n + l,rt) = a horizontal step followed by a single descent 

or a single descent followed by a horizontal step 
= two successive descent steps. 
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Note that the structure of the path classes does not actually depend upon n. This is 
a reflection of the underlying spatial homogeneity of these equations. Thus, for the 
purpose of describing the path classes, one can translate n to 0. 



Now applying (2.13) the difference string equations become 



1 

n 



(bl+i-bl) +3< (a n+1 b 2 l+1 + a n b 2 n+1 - a n b 2 n -a n -ib 2 n ) 



= (a n+1 -a n ) (b 2 l+1 +3t(a n+1 +a n )b 2 l+1 ) +3i (b 2 n+2 b 2 l+1 -b 2 n+l b 2 n ) 

where, for this example, we have set the parameter x equal to 1. The coefficient 

is non- vanishing by (2.11 1 and the fact that the partition functions are non- vanishing 



for n sufficiently large. Hence we may divide it out of the second equation to arrive 
at the slightly simpler system 

1 = (bl+i-bl) +3t(bl +1 (a n+ i + a n )-b 2 l (a n + a n -i)) 



n 

= (a n+ i 



)(l+3t(a n+1 +a n ))+3t(b: 



n+2 ' 



2.2. Motzkin Representation of the Toda equations 



We now pass to a more explicit form of of the Toda equations (2.5) in the case 
i = 2i/ + l: 

1 da„ 



N dt 2l ,+i 
1 dbl 



C 



2i/+l\ 



' n+l,n 



c 



2v+l\ 



■ (a n -a n ^)(C 2v+1 ) 1+ (£ 2 - +1 ) +1 ,-(£ 2 » +1 ) 

v y V / n,7i— 1 V /n+l,n— 1 V / n,n- 



(2.17) 



(2.18) 



N dt 2v+1 

To describe this in more detail we will once again specialize to the trivalent case 
(y = X). There are two relevant path classes here: 

■p 3 (n + l,n)described in Figures [2~T] and [2~2| for the case n = 
V 3 (n + 2, n)described in Figure [2~3| for the case n = 0. 

The latter case corresponds to what was used in [TD] but for Dyck paths (Motzkin 
paths without any horizontal steps) of length 2v. 







Fig. 2.1. Elements ofT 1 (1,0) with two horizontal steps 



Applying ( |2.13 ), the trivalent Toda equations become 
1 da n 
n dt 

+ + (^+2^+1-^+1^) + (b 2 n+1 b 2 l+1 -b 2 n b 2 n ) 



i a l+ibi+i -a 2 n b 2 n ) + (a n+1 a n b 2 n+l -a n a n -ib 2 n ) + (a 2 n b 2 l+1 - a^_ x 6^) 

(,2 l2 



.1*1 

n dt 



- (a„ - a n - 1 ) [a 2 n b 2 n + a„a„_ i b 2 n +0^^ + b 2 n b 2 n _ 1 + b 2 l+1 b 2 n + b 2 n b 2 n ] 

+ (a n+1 bl +1 b 2 n -a n blb 2 n _ 1 ) + (a n b 2 n+1 b 2 l ~a n ^ 1 b 2 l b 2 l _ 1 ) 
+ (a n -ib n+1 b 2 l -a n - 2 blb 2 n _ 1 ) 
where we have again used the relation x= and then set the parameter x = l. 
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Fig. 2.2. Elements o/P 3 (l,0) with no horizontal steps (Dyck paths) 




Fig. 2.3. Elements of V 3 (2,0) 



3. Continuum Limits A number of discrete variables will appear in the fol- 
lowing discussion as we prepare to make the transition to the continuum limit. Not 
all of these discrete variables will be directly involved in the description of that limit. 
However, in order to avoid confusion, it is perhaps best that we start by briefly de- 
scribing all of these variables and their interrelations as well as their relations with 
other variables and parameters. As indicated at the outset, the positive parameter N 
sets the scale for the potential in the random matrix partition function and, through- 
out this paper, it is taken to be large. The discrete variable n labels the lattice position 
on Z-° that marks, for instance, the n th orthogonal polynomial and recurrence coef- 
ficients, diagonal or sub-diagonal. We also always take n to be large and in fact to be 
of the same order as N. As stated in the Introduction, it is to be understood that as 
n and N tend to oo, they do so in such a way that their ratio 



n 
N 



(3.1) 



remains fixed at a value close to 1. In fact within all subsequent proofs and derivations 
x itself will be held fixed. 

In addition to the global or absolute lattice variable n, we also introduce a local 
or relative lattice variable which we will denote by k. It varies over integers but 
will always be taken to be small in comparison to n and independent of n. We will 
frequently study expressions involving n + k which we will think of as small discrete 
variations around a large value of n. We have already encountered this in the form of 
the difference string or Toda equations. The spatial homogeneity of those equations 
manifests itself in their all having the same form, independent of what n is, while k 
in those equations varies over {— v— 1,..., — l,0,l,...,z/ + l}, as explicitly displayed for 
the trivalent case (v=l) in subsections 2.1 and 2.2 Indeed in what follows it will 
suffice to take v+ 1 << n in order to insure the necessary separation of scales between 
k and n. 



14 



CONT. LIMIT OF TODA LATTICES FOR RANDOM MATRICES 



We introduce a number of other scalings of variables that we will be using. 

si=ar^ii (3.2) 

or more generally for v £ Z + , 

S2v+i=x v ~h 2 »+i (3.3) 
k 

w = l + -. (3.4) 

n 

The variable w is introduced for two reasons. First, to help make some of the subse- 
quent expressions less cumbersome. Which value of k is intended will be clear from 
the context or it will be made explicit. The second reason for introducing this vari- 
able is that it provides the transition to the continuum equations. As we shall see, 
at a certain point in the subsequent arguments w will appear within the arguments 
of coefficient functions of the large n asymptotic expansions. Since these coefficient 
functions are in fact analytic in their arguments, we will take advantage of this fact 
to regard these functions as analytic functions of w regarded as a continuous variable. 



In Theorem 3.3 we will present a more precise statement and extension of prior 



results on the free energy expansion (1.17) as they will relate to what we do in the re- 
mainder of this paper. However, before getting to that we need to recall a preliminary 
prior result: 

Proposition 3.1. fTZy 

2 2 ((n + ky 1/2 /n+fc\ 1,-1/2 \ 

T n+k,N V-l^v+l) —T n+kn+k I I — I tl, I — I t 2v +l I 



' n+k,n-\-k 



' k\' 1/2 ( kY~ 1/2 ' 
i + -J s 1 ,h + -J s 2u+1 



2 / 1/2 -i^-1/2 \ 



Proof. The first equality follows from an appropriate change of variables in the 
Szego representation (2.8 2.9), the second just applies the definitions (3.2) and (3.3) 
and the third applies (3.4). The change of variables in the Szego representation 
amounts to introducing the re-scaling Xj = -\fx\j , from which we then have 



<7( n ) {+ i . 



J ■ ■ ■ |exp (- r>{ (\^+h„ + ix»- 1/2 \r +1 



= x n2 ' 2 Z^{s uS2u+1 ), 



(3.5) 



where V(A) = JXj<£ W — -M ■ Shifting from n to n + k this becomes 



y(n+k) 



{t\,t2v-\ 



Jn+kf/2 y(n+k) 
^n+k 



-1/2 



Si, 1+- 



t/ — 1/2 



S2i/+1 



The proposition follows immediately from this. □ 
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Remark 3.2. Starting in subsection \S . 1\ we will in general be setting x=l. That is 
because the main focus in this paper is on the structure of the partition function coef- 
ficients as generating functions for map enumeration. However, for (possible, future) 
applications to the statistics of random matrix eigenvalues, the ability to asymptot- 
ically "detune" the matrix size away from the scale of the potential is important. 
Therefore, we have chosen to keep this parameter free up through this preliminary 
subsection. 



(3.6) 



We introduce one more notational definition: 
A^ogT* in (s 1 ,s 2l H-i)= 1 °g^+/^^ 
Theorem 3.3. JMIMW 

i 2 / \ 2 ~2 / - — 1/2 -v-1/2 \ i 

\ogT n+kn+k {si,s 2 v+i)=n w e (w 'sx,w 1 s 2v+ i)+ 

+ e 1 (w- 1 ' 2 s 1 ,w l '- x ' 2 s 2v+l ) + ... 

+ J^w 2 - 2 ^^- 1 ' 2 *^- 1 '^^) + ... (3.7) 

is an asymptotic expansion in n~ 2 , uniformly valid for (s\,S2u+i) €/C= any compact 
subset of (—6,6) x [O,,?^' 9 ') and \k\ < v-\-l with < e for e and 6 sufficiently small. 
(Here si"' 9 ^ is a fixed positive constant depending only on v and g.) Explicitly, this 
means that for each g there is a constant, C g , depending only on v and K, such that 



l ogT 2 +kn+k (s 1 ,s 2v+ i)-n 2 w 2 eo(w 1/2 si,w I/ 1,2 s 2v +i) - . . . 

1 ~2-2<j I 1/2 -v—XII \ 

v a e g (w ' si,w 1 s 2y+1 ) 



7 2 ff -2 



<9f 
- ripg 



for all (si,S2i/+i) € K. and \k\ <v + l. 

a) This expansion may be differentiated term by term in S\,S2u+i with the 
same type of uniformity except that the constant C g will now also de- 
pend on the multi-index of the derivatives. Moreover, the coefficient 
e g (w^ 1 / 2 sijU;" -1 ' 2 s 2u+ \) and its mixed derivatives in (w,s\,s 2u +\), after be- 
ing evaluated at (w,si) = (1,0), are all complex analytic in a disc of radius 
s 1 ^' 9 ^ centered at in the complex S2v+i plane. (Here we introduce w as a 
continuous complex variable replacing w; hence, one may differentiate e g with 
respect to it.) One expects this radius of convergence to be independent of g 
as is known to be true in the case of even weights p]/. 

b) One also has an asymptotic expansion for differences A^logT 2 • 



n.n ' 



A -fe 1 °gTn,n( s l> s 2«+l)=X)^2^2 ( 3 ' 8 ) 
5=0 

El & 2-2q I -1/2 I/-1/2 \ I f k^ 3 

j= i J - V ' 

where we regard w as a continuous complex variable. Once again this expan- 
sion is uniformly valid for (si,S2v+i) G K. and \k\<v + l, by which we mean 
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that for each g there is a constant, D g , depending only on v and K, such that 



9 1 

A fc logr^„( Sl , S2l/+1 )-^ — 



m-2 



2( ff -m)+l 



E 

3=1 



1 &> 



m— 



w 2 - 2m e. t 



(V 1 



(3.9) 
A, 

,2 S • 



< 



77ns expansion may be differentiated term by term in s\,S2v+i preserving 
uniformity. 

Proof. The basic result is that of [5] extended, in [3J, to the case of weights with 
odd dominant power (see ( 1.17[ ) [j]- From this it follows that one has constants 
C g depending only on j/,/C such that 



1 °g < rn+fc,n+k( s i> s 2 I /+i)-(rc+fc) e (w 1 si,w" 1 s 2 „+i)-... 

■ • ■ - ( n + fc)2g -2 e g(^ 1/25 i>' / ' 1/2 ^+i) 



We then rewrite the above equation using n + fc = mD, 



logT 2 +kn+k (s 1 ,s 2 v + i)-n 2 w 2 e (w 1/2 s 1 ,w" 1/2 s 2l ,+i) 



< 



C a 



(n+k) 



2<) 



,2 S -2 



~2-2o / ~.-l/2 ~.v— 1/2 \ 



< 



g 9 



The desired estimate is realized by taking C g — ( . 1 _ e ^ 2g ^_ . 

The ensuing statements of the theorem, in (a), also follow directly from these 
prior results. In particular, the analyticity of e g in its arguments follows from (1.171 
|uj mixed derivatives then yield linear combinations of derivatives of e g with respect 
to its arguments whose coefficients are polynomials in Si,s 2v +i and fractional powers 
of w (which is bounded away from zero). Evaluating at (iy,si) = (1,0) then yields a 
linear combination of derivatives of e g with coefficients that are polynomial is s^u+i- 
By (JiiJ) this linear combination is analytic in a disc as stated in the theorem. 

For (b) observe that a straightforward estimate of the difference of the asymptotic 
expansions for logr^ +fc +fe and logr 2 „ yields 



1 og^+fe,„+fc(si,s 2l , + i)-logr^ n (si,S2^+i)-n 2 w 2 eo(w 1/2 s l ,w l/ 1/2 s 2v +i) 
+e (si,s 2v +i)- ■■■ 



1 



n 2g-2 

We rewrite this as 
Afclogr 2 (si,s 2v +i) 



2—2q l —1/2 v—l 2 \ i / ' 

w y e g (w 1 si,w 1 s 2 u+i) + e g {si,s 2 u+i / 



< 



w=l+± 



2Cg 

,2 S 



g 1 

E -^=2 (^ 2m ^(w- 1 / 2 s 1 ,w^ 1 / 2 s 2u+1 )~e g ( Sl ,s 2u+1 )) 



m=0 



< 



2C n 



? 2 S ' 
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which is 

AfelogT^ n (si,s 2l ,+i) 

-E^f( i +-r 2m ^((i+-)^ 1/2s i I ( i +-r 1/2 ^ + i) 

± — t n £m a i n n n 

-e g (si,s 2 u+i) 



m—0 



< 



2Q, 
7,29 



One now Taylor expands the e m terms centered at w = 1 and evaluated at 1 + £ . For 
notational convenience set 

F m (w, si , s 2v +i) = w 2 ~ 2m e m (w~ 1/2 s 1 , w u ~ 1/2 s 2y+1 ) ■ 

Then this expansion has the form 

2(g-m) + l 



A fe logT^„( Sl , S2iy+1 )-^ — J2 



m—0 



3=1 



1 d 1 
j\ dm 3 



F m (w,sx,s 2v+1 )\ w=1 



+R^(w,s u s 2u+1 )('^j 



2(fl-TO+l) 



< 



2Cg 

,.2g 



where R( m ' (w,s\,s 2l/ +i) denotes the remainder term, of order 2(gr — m + 1) in w, for 
fm- By elementary inequalities one then has 



AfelogT^ n (si,s 2l ,+i) 

5 2( S -m) + l . 

E i \ "* 10 „„2-2m, 

„2ro— 2 



m=0 



,;=i 



j! 9«P 



< 



2C^ + ^ 

,2s „2g 



£|iJ (m) (u>,«l,*2H-l) \k\ 



2(g-m+l) 



m—0 



By Cauchy's remainder theorem for analytic functions one has 



R (m) (w, Sl ,s 2l/+1 ) 



<2 2(g-m+l) . 



d m = max max F m (w,Si,s 2v+1 ) . 

(fli,«a„ + i)e/C|n;-l|=l/2 

(We assume here, as we may, that e<l/2.) Thus statement (b) is established by 
taking 



D g = 2C g +J2d m \2v + 2\ 



2( fl -m+l) 



m=0 



Remark 3.4. We mention here that the variables Sj as defined above differ slightly 
from their usage in related works \1(K [ffi where Sj = —Cjtj for appropriate constants 
Cj>0. Also for comparison with J3j], 53 = — u where u is the weight parameter in 
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that work. We further observe that because of the combinatorial interpretation of 
these generating functions, one can show that e g (s2 V +i) is even in s 2 „+i. Thus odd 
derivatives of e g are odd functions and even derivatives are even. 

Remark 3.5. By collecting terms in (3.9) of the same order in n^ 1 one sees that 
the asymptotic series represented by (3.8) does indeed have a uniformly valid, well 
ordered expansion in inverse powers of n. The precise form of the coefficients in this 
re-summed expansion is not prima facie obvious; however, the results in the remainder 
of this paper will show precisely how these coefficients and those of other similarly 
derived asymptotic series can in fact be determined. A key point for this process may 
already be observed in \3.9\ ; namely, all terms in the re-summation will be in the 
form of differential expressions in the continuous variable w which are then uniformly 
evaluated at w = 1 . In the relevant settings, the coefficients of the inverse powers of n 
may be regarded as hierarchies of differential eguations in w which are to be solved and 
whose solutions are then evaluated at w = 1 in order to yield explicit expressions for 
generating functions (in S2v+i) and similar functions of combinatorial or statistical 
interest. In the remainder of this section we build upon these ideas to derive the 
general form of the various hierarchies of differential equations. In sections^ and^ 
we carry out this process in complete detail for the trivalent case (v = l). To get an 
idea of how the whole strategy comes together the reader might find it useful to browse 
these last two sections before proceeding systematically through the general derivations 
which begin in subsection \3.1\ 

We will make essential use of the Hirota formulas for the Toda variables in their 
original scaling. 
Lemma 3.6. (Hirota) 



a n ,N : 



log 



h 2 

u n,N - 



_}_d_ 
~Ndh 

1 d 2 



! n+l,N 



'n,N 



i d 2 



l__d_ 
N&h 
1 



log 



Z 



(n+l) 
N 



it\fi v 



+1) 



(3.10) 



k'A^^^W), (3-11) 
where the factors of 1 /N are consistent with the energy scaling chosen in the definition 



N 2 dtj ^ 6 N 2 " N 



of [i (l.D 



Proof. From (2.17) and (2.18) one deduces that the Toda equations for v-- 
1 



= are 



da n 



N 



N 
1 



dti 



= b 



n+l,N ' 



h 2 
°n.N 



N dh 



= b n ,N( a n,N- 



■CLn-l,N 



Substituting the fundamental identity (2.11) into the second of these equations one 
has 



1 d , 2 



- 21 0gT^jv-l0gT^_ 



1,N) 



From the Szego representation one has that the t 2 n are simultaneously analytic in 
(ti,t2v+i)- By continuation of (ti,<2i/+i) back to (0,0), the recurrence coefficients 
become those of the Hermite polynomials. From these initial values and the previous 
line one may deduce that, in fact, 

, .„ r n+l.N 



1 d 



N dh 



'n,JV 
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which is the first Hirota relation. Substituting this into the first Toda equation above 
one may similarly derive the second Hirota relation. The expression in terms of the 
partition function in each case follows directly from (1.16 1. □ 
Corollary 3.7. 



1 d 



a n +k,N(ti,hv+i) — -Jjwr [l°g T n+fe+i,iv(*i)*2i/+i)-logr n+fci j V (ti,t2i/+i)] 



x 1 ' 2 d 
i d 2 



Ailogr, 



n+k,n+k 



(si,S 2v 



+1) 



(3.12) 



&n+k,iv(*l>*2i/+l) — Jp~Qfi \°g T n+k,N 



ri2~Q^2 [ 1 °g r n,n(si> s 2 l ,+l)+A fc logT^„(si,S2„ + l)] . (3.13) 



Proof. These representations follow directly from the Hirota relations (3.101 



(3.111, along with (3.2|, (3.3) and the definition (3.6) of A k . □ 



We have the following asymptotic expansions for a n +k,N, b n +k,N- Moreover, 
expanding the d/dw and d/dsi derivatives within the coefficients in these asymptotic 
expansions one can see that these coefficients acquire a self-similar scaling. 
Theorem 3.8. The following are asymptotic series in 1/n: 



(3.14) 



a n +k,N = h(s 1 ,s 2 „+x,w) = x 1/2 ^2h g (s 1 ,s 2v +i,w) 

g>0 

h g (s 1 ,s 2l ; + i,w) = -w 1 ^ 9 x 

v , 1 & +1 

^— J j\ dsidwi 
2.9i +j = .9 + 1' 
.9i>0,j>0 



,2-2g 



^e gi {{ww)- 1 / 2 s 1 ,{wwY- 1 ' 2 s 2v+l ) 



(3.15) 



b l+k,N =f{sl,S2u+l,w)=X^2f g (s 1) S 2 u+l,U>)n 29 

g>0 

d 2 

fg{si,S2v+l,w)=W 2 ~ 29 ^e g (w- 1/2 S ll W 1/ ~ 1/2 S 2 v+l)- 

Moreover 



hg(si,S2 V +i,w)=w* 9 u g (s 1 w 1/2 ,s 2 v+iw l/ 2) 



fg(si,S 2v+ i,w) ^W 1 29 Z g (s X W 1/ ' 2 ,S 2u+ iW U a) 



1/2 



(3.16) 
(3.17) 



where u g and z g are analytic functions of their arguments in a neighborhood of (0,0) 
and w is a continuous variable in terms of which the general form of these coefficient 
functions is described.. 



Proof. We first consider (3.15). By (3.13), (3.7) and the fact that these asymptotic 
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series may be differentiated term by term, one has 

d 2 



h 2 

°n+k,N ■ 



1 a 



^— ' n 2 9 ds 2 

3=0 1 
d 2 



9=0 



Thus we dehne 



d 2 

z g {siw~ 1/2 iS 2 v + iw v ~i) = -g-^e g (q 1 w u ~ 1/2 s 2 v + i)\ q = w -i/2 Sl 



v+l) |g=w _1 / 2 si • 



(3.18) 



which, in light of Theorem 3.3 establishes all claims concerning / and f g . The case 
for h and h g proceeds in essentially the same manner but is a bit more complicated. 
By (3.12) and (3.7) one has 

x i/2 Q 



&n+k,N : 



n ds 



- ]T [(n + k + l) 2 - 2 % ((w+ l/n)-^ 2 Sl , (w + l/nr- 1/2 .s 2 „ +1 ) 

1 S>0 



in + kf-^e g (w-^s 1 ,w v - 1 /h 2u+1 ) 



Setting w = 1 + this may be rewritten as 

x 1 / 2 d 



h{s 1 ,s 2v +i,w)- 



n ds\ 



9>0 



J2(n + k) 2 - 2 ° w 2 - 2 °e g ((wwr 1 / 2 s 1 ,(wwr' 1 / 2 S ^ +1 



-e g \w 1/2 s 1 ,w 1 ' 1/2 s 2 u+i 



We next expand the summands in terms of Taylor series in the continuous variable w 
centered at w = 1 and evaluated at w = 1 H — ]~r ■ 

n+k 

aji/2 q 

h(s 1 ,s 2u+ i,w) = — — 22(n + k) 2 - 29 x 



n dsi 
1 di 



3>0 



^{n + k) 1 -^ 



3>1 

n dsi 



w =i (n + ky 



g>0 



E 



2^1+ J =.9 + 1 
.9i>0,j>0 



1 d j 
j\ dvji 



w 2 - 2 ^e gi [[ww)- 1 l 2 s l AwwY- 1 / 2 S2v+l ) 



9>0 



E 



2<Ji +i =.9 + 1 
.9i>0,j>0 



i 

j! dsidwi 



„2-2g 



1 e 9l ((^)- 1 / 2 Sl ,(^)^ 1 / 2 S2y+1 ) 
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In the second equality above, we have collected the coefficients of (n + fc) 1-5 ; and 
in the last equality, we make use of the relation n + k = nw. To see the self-similar 
structure of the internal sum in this last line, observe that 



d_ 
dw 



+ ( 1 /^K- 1 (H V - 1 ^ 2 , + l)^^)- 1/2 Sl,«2) 92=M -'/^ +1 



where E(qi,q 2 ) is an arbitrary analytic function of (51,52)- We see from this equation 
that a w derivative of a power of w times a function of the self-similar variables qi = 
(ot)' 1,,2 si, q 2 = {ww) l/ ~ 1 / 2 S2v+\ has the same form, with the power of the pre-factor 
reduced by 1. Thus, by induction, the summands of the expression for h g (si,s 2u+ i,w) 
are of the form 



dsidwi 



W — 1 



„3/2-2 S 



(ww) 1 / 2 



w—1 



= wV*-*E j (w-V 2 8 1 ,w''- 1 '*s 2v+1 ), 



for a function Ej of the self-similar variables w 1 / 2 Si and w v x l 2 s 2v +\. The claims 
concerning h and h g now follow from these observations. □ 

Example. The terms of order less than 1/n 2 in these series have the coefficients 



Q2 

h {si,s 2v+ i,w) = -w - w 2 e {{ww)- 1/2 s 1 ,(wwy- 1/2 s 2v+1 ) 
osiow 



d 2 



ds\d{ww) 

d 2 



-(wwfeoiiwwr^s^wwy-^s^) 



dsidw 



w 2 e (w- 1 l 2 Sl ,w^ 1 / 2 s 2v+1 ) 



(3.19) 



at order 1 in h, where in the second line we have made a change of variables while in 
the third line we have relabeled ww as w: and, 



d 2 

fo(si,s 2ly+1 ,w)=w 2 — ^e (w~ 1/2 s 1 ,w u ~ 1/2 s 2 ,y + i) 
os 1 

1 d 3 

h 1 (si,s 2l/+1 ,w) 



2 dsidw 2 
1 d 3 



w 2 e {{ww) 1 / 2 si,(wm;) 



w—l 



2 dsid(ww)' 



;(ww) 2 e ((ww)- 1/2 Sl ,(wwy-^ 2 ) 



1 d 3 



V2 S1 ,,,— 1/2^ 



~eo(w ' si,w 



2 dsidw 2 

at order 1 in / and order 1/n in h respectively 



(3.20) 
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In particular, the terms up through order 1/n of a nj jv and b n ,N are, respectively, 
x 1/2 h (s 1 ,S2 U+1 ,l)^x 1/2 u (s 1 ,s 2u+1 ) 

= -^ 1/2 7T^r^ 2 eo(^- 1/2 si,^- 1/2 S2,+i)U=i (3.21) 
asiow 



x 



1/2 h 1 (s 1: S2v+i 1 l)=x 1/2 u 1 (s ll s 2v +i) 



1 r) 3 

U ,2 ^- 2 w 2 e,{w-V 2 Sl , W ^l 2 s 2 „ +1 )\ w=1 ( 3. 22) 



■-x-^e (si,s 2u+1 ). (3.23) 



z/o (si , s 2v +i , 1) = £2o (si , s 2 „+i ) 
ds 2 ' 

In Subsection |4 . 1 . 1 1 we will show that the coefficients Uo and z defined here are indeed 
the same as the functions introduced in Section [T] to describe how the endpoints of 
the support of the equilibrium measure depend on the parameters in the exponential 
weight. 

We also introduce a shorthand notation to denote the expansion of the coefficients 
of h(si,S2u+i,u>) and f(sx,S2 V +i,w) around w = l. This is analogous to what was done 
in interpreting (3.8) via (3.9), the main difference being that the order of summation 
is interchanged. This is again justified by the asymptotic interpretation (3.9) where 
both summations are finite. 

Definition 3.9. For w = l + k/n with \k\ <2v and — <e, 

h.,,(m) L.=i / k 



* — ' ml \ n I 

m=0 v 7 

f( Sl ,s 2v+1 ,w)= f; (-Y (3.25) 



mi 

m— u 

where the subscript u/™ 1 -* denotes the formal operation of taking the m th derivative 
with respect to w of each coefficient of h (respectively f): 

^ gm ^ 

h w M = g^ h g ( Sl , 82v +i, w) — 

g>0 

gm ^ 

= E g~^ h (*i , ,w)-^ 

As valid asymptotic expansions these representations denote the asymptotic series 
whose successive terms are gotten by collecting all terms with a common power of 1/n 
in {3.2$ (respectively 

In what follows, in the rest of section 3 and in section 4, we will frequently abuse 
notation and drop the evaluation at w = l. In particular, with x=l, we will write 

m=0 x ' m=0 g>0 



& «+^=E^f (~) =E^,E^/ 9 ( s i'^+i^)^(^) (3.27) 

m=0 v 7 m=0 g>0 











1 




n 2g 
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In doing this these series must now be regarded as formal but whose orders are still 
defined by collecting all terms with a common power of 1/n. They will be substituted 
into the difference string and the Toda equations to derive the respective continuum 
equations. At any point in this process, if one evaluates these expressions at w = 1 
one may recover valid asymptotic expansions in which the a n +k,N and b^ +k N have 
their original significance as valid asymptotic expansions of the recursion coefficients. 
In particular, in Section 5, the results of the formal derivations will be evaluated at 
w = 1 and we will recover explicit expressions for the e g appearing in the asymptotic 
expansion of the partition function. 



3.1. Continuum Limits of the Difference String Equations We are now 

in a position to substitute our asymptotic expansions for a^k^ arid b n +k,N into 
the difference string equations, (2.15) and (2.16). Collecting terms in these equations 
order by order in powers of 1/n we will have a hierarchy of equations that, in principle, 
allows one to recursively determine the coefficients of (3.14) and (3.15). We will refer 
to this hierarchy as the Continuum Difference String Equations. (Note that one has 
such a hierarchy for each value of v.) Of course this is a standard procedure in 
perturbation theory. The equations we will derive are odes in which w, now regarded 
as a continuous variable, is the independent variable. The variables s\ and S2 V +i here 
are parameters on which the ode depends analytically. One must still determine, at 
each level of the hierarchy, which solution of the ode is the one that corresponds to the 
expressions given for h g and f g in Theorem 3.8 This amounts to a kind of solvability 



condition which will be imposed through a small number of initial Taylor coefficients 
of e g to insure that the solution coincides with its enumerative interpretation in terms 
of counting maps. This solvability analysis will be illustrated in detail in Section [4] In 
this subsection the main emphasis will be to derive the form of the continuum string 
difference equations and their general solutions. 



From this point on in this section (and in fact for the remainder of the paper) 
we will set x=l; i.e., n = N. This has the effect of centering the matrix size n at 
the same scale as that of the potential, N. We will also set s± — from now on 
since its role in determining the structure of the asymptotic expansions of a n +fe and 
b n +k is now completed. When x=l, S2v+x =~t2v+x] however, we will continue to 
present statements in terms of the s-variables. If one wants to subsequently "detune" 
to a value x<l one can do this by replacing s 2 „+i with its expressions in (3.3) 
and comparing to (3.14) and (3.15). Finally, when the context is clear, we will for 
simplicity just use s to denote S2u+i- 



We begin by substituting the expansions (3.26) and (3.27) into the difference 
equations, (2.15) and (2.16), satisfied by these coefficients (as represented through 
C). We arrive at the following formal asymptotic equations. For equation (2.15) one 
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has: 



E 



h w(m) ( 1 



+ (2v + l)s 2 u+l 
2/i(P)+l oo 



E ; 

_m— 



E 

PG-P 2 "(1,0) 



^ to! I n 

p a — 1 m— 



n e^ 

Pb — 1 m— 



iH / «Pb 



m! \ n 



E 

PG-P 2 "(2,0) 



2M(P) oo 

n\ ^ fWm) / lp a 

^ m! V n 

p a — l rn—0 v 



y-M(P) + l oo 

n e 



A, 



2 M (P) 



n e- 

p a — 1 m— 



o('«) / «p, 



4. - 1 



Pb— 1 m= 

i/-/*(P)+l oo 



TO! \ TT, 



C"») / «Pb 



(3.28) 



n e 

Pb — 1 rn— 



/ — 1 



(3.29) 



where /i(P)= [f/2j for <r equal to the total number of horizontal steps in a given 
path P and £ Pa (respectively £ Pb ) denotes the lattice location of the path at the 
p 1 ^ horizontal step (respectively before the pjf 1 downstep). Note also that we have 
taken advantage of the discrete space homogeneity of these walks in order to shift the 
initial/final points of these paths to (1,0), (2,0) and (2,1) in the respective cases. Note 
further that each term of (3.28) is divisible by 



oo 

to! 

m— 



Likewise for (2.161 we find 



1 

n 



E 



fu 



+ (2v + l)s 2v+x 

2/i(P) + l oo 



E 

PG-P 2 "(1,0) 



n e 

p a — l m—0 
2a*(P)+1 .oo 

n E 

p a — 1 m—0 



h w{m) ( I 



TO! \ TT, 



fu 



TO! \ TT 



Km ftp*- 1 



V— M(P) oo 

n e 

Pb = l m—0 

u—fj,(P) oo 

n e 

Pb = l m—0 



(m) I l p 



(3.30) 
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To begin with, the equations at leading order are 



(2v + l)s 2 „ +1 J2 



2v 

2(1,1/ — (1 — 1,1/ — jJL- 



J^(«^ +1 )//o, 



l = d w fo + {2v+l)s 2v+l ^ 



/i=0 



2// 



— fJL — 1,1/ — (i 



u w I JO 



This may be written in a vector form as 
Proposition 3.10. 



ho w 



/o,i 



where 

v-1 

A u = l + (2!/ + l)s 2l , + i^ 



2;/ 

2fJL + l,V — fl—l,U — 



(2iy + l)s 2u+1 J2 



/i=0 



2z/ 

2(1,1/ — jJL— 1,1/ — (1- 



A 12 = (2is+l)s 2v+ i^2 



/i=0 



2// 



2(j,,i/ — [i — l,i/ — (X- 



\ t 2/1+1 fi>-n-\ 

1 )(2/x)^- 1 /o^, 
1 )(^-/x + l)^/o-' 1 - 1 ) 



A-21 — /o-A-12 , 

I/-1 

A 22 = l + (2z/ + l)s 2l/+ i^ 



/i=0 



2i/ 

2(i+l,i/ — (i — l,i/ — (i 



(3.31) 



(3.32) 



(3.33) 

(3.34) 

: Au, 
(3.35) 



in terms of tri-nomial coefficients. 

We will see that the coefficient matrix of the n~ 2s_1 terms is of the same form 



as the matrix in equation (3.31). Thus the following lemma will be useful: 



Lemma 3.11. For the matrix A given in (3.31) 

^-i _ I fo,w ^o, 
ho,wfo fo.i 



Proof. The second column of the inverse follows directly from (3.31 ). To find the 
first column one notes that 



An A 



12 



1 



An -A 



12 



A12/0A11/ A 2 n -A 2 12 f \-Aufo An 

Thus we have that An/det(A) = fo. w and — Ai 2 /det(A) = ho tW , and the result follows. 
□ 
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We will also need the following lemma: 



Lemma 3.12. For the matrix A given in (3.311 

tf ho Ai2-5/ o A u = 
<^ o A 2 2-<i>/ o A 21 =0, 

where Sh (resp. Sf ) denotes the functional derrivative with respect to ho (resp. fo). 
Proof. The left-hand side of the first equation, written out, is 



(2z/ + l)s 2 H 



E 



2v 



' \2n,v — f/, — l,v — n + 

v—l 

E 



li=0 
u-1 

E 



2v 



2fi+l t u — fj,— l,v — 



2v 



1 )(^-M + l)(2^)^" 1 /o" A ' 1 

)(u-^-l)h 2 / +1 fr"- 2 (3-36) 

1,V j-l/ 



= x \2(i,v-[j,-l,v-n- 



2fi— 1 pv-^i—l 



\{2n){v-n)hf^t Q 



Shifting the index of the middle sum by one sees that the coefficient of each 

monomial in h ,f cancels and the result follows. To prove the second formula one 
first applies the identity (3.34) to find 



5 h A 22 - Sf A 2 i = 5 ho A 22 - foSfo A12 - A12 . 



(3.37) 



Written out, the right-hand side of (3.37) is 

(2v + l)s 2 „ +1 J2 L , , n 

^ \_\2^ + l,v — /i-l,z/-/x 

2;y 

2ji,v — n~ — /i + 1 
2^ 

2//,!/ — /X — 1,V — (1 + 1 



(i/-/i)(2/x + l) 



v—\x—\ 



whose coefficients manifestly vanish. □ 

The homogenous terms of the equations at level n _2£,_1 can be computed directly. 
They are linear in ft, 2g and / s with coefficients depending only on ho, fo and their w 
derivatives. The inhomogeneous (forcing) terms depend on hj for j<2g, and fj for 
j < g. As usual in perturbation theory, the homogeneous part of the equation can be 
derived by replacing (h ,fo) in the leading order equations with 

(h + eh 2g ,f + ef g ), (3.38) 

and retaining just the first order in e terms. We find that the homogeneous terms are 



h-2g,w 
fg,w 



<-0,w 



O.il' 



■fg S fo A 



TO.W 



(3.39) 
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We also have the identity 



A. [ h f 

J a 



= A 



S ho A[ < J 

n 0,wjg 



fo,wh2g 
fd.wfg 



(3.40) 



Lemma 3.12 implies that (3.39) is in fact equal to the right-hand side of (3.40). Thus 
the equation at order n _2s_i has the form 



A. [ h f 

J 9 



(3.41) 



where F^ 1 and are expressions invovling the lower order terms in the expansions 
of (|3.28|) and 13.301). We also find, at order n~ 2g , that 



^{AllVll^j! 



(3.42) 



In fact there is a second equation involving h-2 g -i'i however, it must be equivalent 



to the first and so we do not record it. Equations (3.41) and (3.42), together with 



Lemma 3.11 yield the following 

Proposition 3.13. The functions /i2g,^2g-i,/g may be recursively found by the 
formulas 



l 2g _ 



fo,w ho,- 

fg J \ho, w fo fo 

1 

An 



F, 



(1) 



-in 



F, 



(2) \ dw ' 



2y 



l 2g-l 



F^dw. 



(3.43) 
(3.44) 



3.2. The Continuum Limit of the Toda equations Analogous to what 
was done in the previous subsection for the difference string equations, one can study 



the system (2.17) and (2.18) expanded on the formal asymptotic series (3.26) and 
(071: 



1 d 

n ds 



h(s,w)= 

2fJ,(P) oo 



n\r^ h w { m ) ( £ Pa 
2-^ rn\ 



p a — 1 m— 
2p(P) oo 



nsr^ h w { m ) f £ Pa 



p a — 1 m— 



n e^ 

Pb — 1 rn— 



i/-/i(P)+l co 

n e 

Pb — 1 m— 



/ui< m ) ( ^Pb 1 



(3.45) 
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ld_ 

n ds 



f{s,w)= Yl 

Pe{-p 2 "+ 1 (2,o)} 

2(ti(P) + l oo , 



n e- 

p a — 1 m—0 
2 M (P) + 1 oo 

n e- 

p a — 1 m—0 



m! 



i/-/*(P) + l oo 

n e 



fv 



-2 



\m=l 



„(m) 



m! V n 



E 

pe{-p 2 "+ 1 (i,o)} 

-M(P) + 1 oo 

n e 

pt, — 1 m— 



Pb — 1 m—0 
i/-/i(P)+l oo 

n e^ 

Pb — 1 m—0 
2ju(P) oo , 



TO! 



(™) / «P6 



-2 



(3.46) 



n e- 

p a — l m—0 



'ii>< m > / 



„(m) / « p , 



TO! 



where, again, /z(P) = L (J /2J for cr equal to the total number of horizontal steps in a 
given path P and £ Pa (respectively l Vb ) denotes the lattice location of the path at the 
p 1 ^ horizontal step (respectively before the pl h downstep). As before, the expansion 
order by order produces a hierarchy of equations that we will call the Continuum 
Toda equations. At leading order in the hierarchy one has, for general u, 



ds ^ o \2fi + l,v- fj,-l,v- /J. + 1J 



(3.47) 
(3.48) 



At orders n 2g the equations are equivalent to a hierarchy of 2 x 2 quasi-linear 
systems of PDE: 
Proposition 3.14. 



d (h 



ds V /, 



'2(, 



B 



2tf 







-rifo, w h2g + r 2 hQ^ w f g J + yForcing[ 



( Forcing[ 



where 



(3.49) 



B n — b 22 =y~^ 



2i/ + l 



1 )2M^- 1 /o" 



B-t( 2( ,,„_ 2 ::^ +1 )<-^'>c/.-«, 



(3.50) 
(3.51) 
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B21 — /c3i2 , 



v-1 



r-i - 



2u + l 



2v+l 

2[i + l, v— (i — \,v — jtt+1 



^i\2fi,v-fi,u-fi + lj 



(3.52) 

(3.53) 
(3.54) 



Proof. The proposition follows from the same approach and methods used to 
establish Proposition |3.13| The various relationships between the entries of B are 
consequences of simple trinomial identities as was the case for the relations between 
entries of A. □ 



Remark 3.15. We observe that the homogeneous terms in (3.4-9) are almost in pure 
conservation law form, with the first terms on the right being in the form of a spatial 
derivative of a flux pair associated to the density pair (h,2g,f g )- The form of the non- 
conservative homogeneous terms (the second group of terms) suggests that one might 
be able to use the difference string equations to rewrite these in terms of lower genus 
expressions and thereby pass them into the forcing so that the equations would then 
have the structure of a hierarchy of forced conservation laws. Indeed, in Section^ we 
show that this is what happens for the genus 1 equations. 

We also note that the numerical coefficients appearing in the homogeneous terms 



of (3.4-9) depend only on the total number of Motzkin paths of each class V J (rn\,m,2) 
that appear in the Toda equations. On the other hand the specific terms in the forcing 
expressions depend on more detailed combinatorial characteristics of these Motzkin 
paths. See \1(A 0? for details about how the forcing terms can be determined explicitly 
from the structure of the lattice paths in the case of even valence. 



4. Specialization to the Trivalent Case 

We now illustrate the results of section [3j in the trivalent case (when v=l) to 
demonstrate their form and utility. We will then pick up from lemma |3.13| in this 
case and find explicit expressions for h\,h<i and f\ in terms of /o, h§ and their w- 
derivatives. As mentioned at the start of subsection [3~T| this will require comparison 
with the enumerative significance of the coefficients in the asymptotic expansions in 
order to determine a unique solution. Note in particular that we will need to go 
beyond the results of section [3] in order to explicitly determine the forcing terms F^\ 
and . In this section we will use s to denote S3. 

4.1. String difference Equations for j = 3 In the trivalent case one has 

= bl +1 [(a n+1 -a n )(l + 3s(a n + a n+1 ))+3s(bl +2 -bl)] , (4.1) 
~ = K +1 +isbl +1 (a n + a n+l )) - (bl + Zsbl(a n +a n ^)) (4.2) 
= (b 2 n+1 - b 2 n )+3sbl +1 (a n + a n+1 )- 3sb 2 n (a n +a n _i) . (4.3) 



Using (3.10) and (3.11), and Taylor expansions around w = l of the continuum 
limits we have: 
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Dividing {11]) by b 2 n+1 : 

= (h(s,l + -)-h(s,l))(l + 3s{h(s,l) + h(s,l + -))+3s(f(s,l + -)-f(s,l)) 



E 



?; 

h w (m) 1 



in 



n" 



l + 3s 2h- 



E 




where the second line is evaluated at w = 1. 
• The equation (4.2) becomes: 

- = (/(*, 1 + -) - f(s, 1)) + 3s/(s, 1 + 1) + 1 + -)) 



-3a/(«,l)(/i(a,l) + fe(s,l--)) 

n 




where the second line is evaluated at w = 1. 



4.1.1. Leading Order 

The leading order of the system ( 4.4|4.5 1 is as given in (3.31 ) with v = 1: 



1 + 6s/io 
6s/ 1 



6s 

-6s/io 



(4.6) 



Expanding we have 



= (1 + 6s/i )/io,tu + 6sf 0iW 

1 = 6sf h ,w + (1 + 6s/i )/o 
These can be anti-differentiated with respect to w: 

Ci(s)=ho 
w + C 2 (s) = 6sf h + f 



3shl + 6sf 



The Ci(s) and C%(s) are constants of integration which must be determined by the 
combinatorial interpretation or some other constraints. For example converting to the 
self-similar variables, s = w 1 ^ 2 s, dividing the first equation by w 1 ^ 2 and the second by 
to, we find 

w^ 1 ^ 2 C\(s) = uo(s) + 3suq(s) 2 + 6szo(s) 
1 + w~ 1 C 2 (s) ~6sz (s)u (s) + z (s). 

In order for the left hand side of these equations to give functions of the self-similar 
variable s we must have that Ci(s) = c\s~ 1 and C2(s) — c 2 s^ 2 . However the functions 
on the right hand side are analytic in a neighborhood of s — S3 = and so we conclude 
that Ci—c 2 = 0. Comparing these with (1.10 1 and ( |1.11| , we have shown here that 
the leading order functions of the asymptotic expansion of a n ^N and & ni jv agree with 
the functions uq and zq describing the equilibrium measure. 
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4.1.2. rT 2a terms 

The odd terms of the expansion for h(s,w) are governed by either of the n~ 2g 
terms of equations ( 4.4|4.5 1; i.e., by either of the equations 

= h 2g -i iW [1 + 6sh ] + h 0tW [6sh 2g -i\ - F^_ x 



(2) 

= 6sf h 2g -i iW +Gsfo, w h 2g -i - F^ g '_i , 



4.5) which 
cnt to the 



where the f 2 g-i are the forcing terms coming from the n~ 2g terms of (4.4 
do not contain a n ^23 -1 or its derrivatives. The first equation is equiva 
specialization of (3.42) to the trivalent case: 

d w {[l + 6sh Q ]h 2g _ 1 } = F 2 {1 g ) _ 1 . 

4.1.3. n 25 1 terms 

The even terms of the expansion for h(s,w) and the terms of the expansion for 
f(s,w) are governed by the n _2fl_1 terms of equations ( 4.4||4.5 ). We find the system 
(given by taking v — 1 in (3.41 )) 



l + 6sh Q 6s \fh 2g \ l \_ 
6s f l + 6sh I /, 



E 



F, 



2y 
'2) 



(4.7) 



2y 



where F 2 ) are the forcing terms coming from the n 29 1 terms of (4.4|4.5 1 which do 



not contain an h 2g or f g or their derrivatives. 

Applying Lemma |3.11| and Proposition |3.13| we have 



>2y 



ho,wfo fo,-< 



F, 



(i) 



F, 



X \dw. 



(2 



2;/ 



l ^- i= ukhJ F2 ^ dw - 



(4.8) 
(4.9) 



Lemma 4.1. hi( s,w) = ^hp tW (s,w) 



Proof. From (3.44) and the n coefficients in the first difference string equation 



(4.4) one has 

[l + 6s/i ]/ii ; 



(l + 6sh )h , w +6sf , w +C(s) 



({l + 6sh )h 0iW ) + C(s)by g^l 
C(s) 



hi — ^.h 

2 1 + bsho 

We see that this agrees with the first two terms of the asymptotic expansion of a n +k,N, 
given in ( [3Tl9| ) and ( [3T20] ), 

h (s,w) 
hi(s,w) = - 



d 2 



'1 Oh j Oil' 2 U ( 



dsidw 
1 d 3 



si=0 



w 2 e Q ( 



1/2 Sl ^ 2 s 3 



(4.10) 
(4.11) 



si=0 
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from which we can also conclude that C(s) = 0. 
Proposition 4.2. 



fo.w h . 

ho, w fo fo,; 







Proof. From (3.43) and the n 3 coefficients in the differenced string equations 
( |4.4||4.5[ ) one has 



fo,w ho,- 
ho, w fo fo,i 

6s(hi + h ^ w ) (h\ tW + \ h ww } + (1 + 6sh ) {^hi,ww + \h 0www } + 

. www ,w ^0, ww 

,www 



dw. 



Substituting for the h\ on the right-hand side, from the previous lemma, one finds 
that the integrand is an exact derivative so that the right-hand side becomes: 



fo,w ho,- 
ho, w fo fo,-, 



\rs(ho, w ) 2 + 2 s hoho,ww + 4s/o iJOTt , + -jqho^ ww + Ci(s) 

2s^o,w/o,«) + s {hofo ,WW ft'O.ww fo) + afo,ww + C 2 (s) 



(4.12) 



One finds further, using the leading order equations ( |4.6[ ), that the second entry of the 
right vector in (4.12| is identically zero, except possibly for the constant of integration 
C 2 (s). 

Converting to equations in the self-similar variable s = w 1 / 2 s, and focusing on the 
terms involving C\ and C 2 in the expression for f\ (the second component in (4.12|), 
we have 

fi(s) = w~ 1 z 1 (s) 

= [terms not involving C 1 ,C 2 ] - {(/o Vw)Ci(s) + (fo, w )C 2 (s)} 

z i(s) — [terms not involving C\,C 2 \ — -zo(s)(uo(s) + su (s))w 3 / 2 Ci(s) 

-(zo(S) + ^Sz , (S))wC 2 (s). 

For this to give an equation for z\ as a function of the self-similar variable s we must 
have that Ci(s) = C3S 3 and C 2 (s) = c 2 s 2 . To pin down ci and c 2 we expand the first 
3 terms of the Taylor series for z±, as given just above, and find that: 



zi(s): 



-c 2 s 



-72c 2 + 6c 3 + 810)s 4 



On the other hand from the asymptotic expansion (3.18) we have 

d 2 



zi(s) = /i(s,l) = g^eiis^s) 



(4.13) 



(4.14) 



s,=0 



and thus the combinatorial meaning of the jth coefficient in the Taylor expansion 
of z\ is the number of genus 1 maps, with 2 vertices of valence 1, and j vertices of 
valence 3. The set of maps with a fixed valence structure on their vertices is bijectively 
equivalent to the set of pairs of permutations (cu,a), where a is a fixed permutation 
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whose cycle structure matches the valence structure of the vertices and uj is a fixed 
point free product of disjoint transpositions, satisfying a further condition equivalent 
to connectedness of the corresponding maps. It is straightforward to partition these 
pairs by the genus. More details on this equivalence can be found in [T] and |16j . 
This allows one to efficiently count the number of maps corresponding to the first few 
Taylor coefficients of z\ . The numbers of genus 1 trivalent maps for j — 2 and 4 are 
found to be and 810 • 4! = 19440 respectively. Therefore we conclude that c 2 = c 3 = 0. 
□ 

4.2. Toda Equations for j = 3 

In the trivalent case we find that the leading order equations become 

-^ s h (s,w) = 3d w [f$ + h 2 fo] (4.15) 

-^fo(s,w) = W w (hof$) +3(d w h Q ) [/ 2 + ^/o] (4.16) 

One could integrate these equations and, after determining the constants of integra- 
tion, show that they are equivalent when w = l to (1.10) and ( |1.11[ ). As this has 



already been done for the leading order of the continuum difference string equations 
in section [4. 1.1 1 we omit the analogous computation here. 
The higher order equations (for g > 0) are 



d 
ds 3 



l2g 



= 39,, 



2h f (2/ 

/ (2/ + /l§) 2h f 



-3(2/0+^ 







howfg 



fowh<2g 



I2g 

h 

(Forcing' 1 - 1 
Forcing^ 2 - 1 



Remark 4.3. We now observe that, for the case of g = l, by using Proposition \4-S\ 
to re-express the terms in the second summand above in terms of h ,f and their 
w -derivatives, these terms may be absorbed into the forcing and those homogeneous 
terms that remain are now in pure conservation law form as was asserted in remark 

mm 

4.2.1. Odd Terms 

The odd terms of the expansion of h(s,w) also generate a hierarchy of (scalar) 
quasi-linear pde, which are recurisevly decoupled from the even terms. The odd terms 
do appear in the forcing terms for the non-homogeneous equations determining /i 2g 
and f g described in the previous subsection. 



The n 29+1 term of the expansion of (3.45) is 



= 30 w (2/ioH-i/o) + Forcing 2g _ 1 . (4.17) 



5. Determining e g 

Recalling the basic identity (2.11) 
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we have, by taking logarithms, 



logr 2 +1 - 21ogT„ 2 + logT„ 2 „ 1 =log(^ ) _log(fe2 )( ) , 



(5.2) 



where the initial value &„(Q) = n is given by the recursion relations of the Hermite 



polynomials. As in [10] , we can use formula (5.2) to recursively determine e g in terms 
of solutions to the continuum equations. We use the asymptotic expansion of b 2 n which 



has the form (3.15): 



1 oo 



-2fl 



3=0 



3=0 



where we have used the self-similar scaling: 



f g (s,w)=w x 29 z g (sw 1/2 ), giving f g (s,l)=z g (s). 



(5.3) 



(5.4) 



In this section, unless otherwise stated we will use s to denote s 3 . It should also 



be noted that the left hand side of equation (5.2) has the form of a centered second 
difference, Air 2 „ - A_ir 2 n . 

We introduce here the classes of iterated integrals of rational functions (or iir for 
short). These classes are defined inductively in terms of the variable z = zq regarded 
as an independent variable. To begin with, the class contains rational functions of 
z. One then adds integrals of these rational functions with respect to dz. Next one 
considers the vector space of polynomials in products of these integrals over the field 
of rational functions of z and augments the space by integrals, with respect to dz, of 
these functions. One continues this iterative process up to any given finite stage. (In 
the classical literature these classes are sometimes referred to as abelian functions.) 

It follows recursively from (3.43) and (3.44) that u 2g ,u 2g ~i and z g are possibly 
in a larger class of functions given by iterated integrals of rational functions as well 
as square roots of rational functions of Zq. 

For general g we have the theorem: 
Theorem 5.1. The function e g (ss) satisfies 



e ff (s 3 ) : 



72 + 1 



,-71—1 



s 71 H g (s)ds — s 3 



— ti— 72— 2 



s 7l+72+1 ff 9 (s)ds 



(5.5) 



where H g (s) is a collection of recursively defined drivers for e g involving terms de- 
pending on Zj for j <g and e, for j < g, with 



(71,72) = (1-4 5) 1) or (3-40,-3), 



(5.6) 



and where C\ and C 2 are constants of integration determined either by the analyticity 
of e g {s^) or the initial Taylor coefficients of e g (s3) determined by some other method 
(for instance direct counting of maps with few vertices). Either choice for the pair 
(71,72) in (5.6) produce the same expression. Moreover, assuming that Zj for j<g 



are class iir, e s (s3) is also in the class of iterated integrals of rational functions. 



Here s is a variable of integration, although it plays, in the integrand, the role of 



S3- 
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Proof. We start from the expression (5.2) and inductively assume that all neces- 
sary Zj have been determined (we will need them for j<g, and will also need that, 
by induction, ej has been determined for j <g). The left hand side of (5.2) is a sec- 
ond order centered difference and so has an expansion for large n involving only even 
derivatives of the spatial variable w. We have 



y— 

g>0 



dw' 



-w 



2-2g 



+ 



i=i 



1 d 4 

Udw 4 
■3 Z J-1 Z 1 



d 2g 



(2g)\dw 2 9 

(-iy +1 z{ 



w 2 e 



dog(z )+ 



'0 



(5.7) 



where on the left hand side eh = eh{w x / 2 s^). 

In the coefficient of n~ 2g on the left hand side one expands the term containing 
e g as a second order linear differential operator applied to e g . The remaining con- 
tributions in the equation are terms that have been recursively determined. More 
precisely, expanding the left hand side of 



dw'- 



we find 



(2-2g)(l-2g)e g - 
We then multiply by S3 1 where 

and integrate once to find 



-(7-8 5 ) S3 e; + - S K = ^( s 3)- 



7i = l-4ff,3-40 



(2 -2.g) (l-2 5 ) +2 



71 + 1 

Next multiply both sides by S3 2 with 



--s 7l+ V • 



s 11 H g {s)ds + C' l . 



72 = 1,-3, 

respectively for each choice of 71, and integrate once to find 



1 



,71+72+2 _ 



7i 



H g (s)dsds' + C 1 sf +1 +C 2 



(5.8) 

(5.9) 
(5.10) 
(5.11) 
(5.12) 

(5.13) 



We conclude by switching the order of integration in the double integral and comput- 
ing the integral with respect to s'. 

Finally we note that the H g (s3) are functions of Zj for j < g, and ej for j < g, and 
so provided that the Zj are in the class of iir functions, we have recursively that H g is 
an iir function. Therefore a consequence of formula (5.5) and this assumption is that 
e g will be an iir function. 



Corollary 5.2. If the driver terms have Taylor expansion 

00 

ff S (s3) = $> S »< 



. 2k 



(5.14) 



k=l 



3G 
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then the Taylor coefficients of the e g (s3) take the form 

.{3),„^ _,o^ ( 2fc ) ! 



«W (2k) = Tin (2k) -, r- 

9 y ' WK ' (l-2g + k){2-2g + k) 



(5.15) 



fork^2g-l,2g-2. 

The proof follows from integration using the power rule. The exceptions, or pos- 
sible resonances, at k = 2g— 1 and 2g — 2 occur at precisely the powers of s associated 
with the constants of integration C\ and C2 ■ To demonstrate the usefulness of Theo- 
rem |5.1[ we will now compute explicity expressions for eo and e\ as functions of the 
fundamental auxiliary variable z$. These are analogous to the expressions determined 
for the e g in [TP] , 

5.1. Example: g = 



In the case of g = 0, we have that Hq(s) = log(zo), and Theorem 5.1 gives 

fS3 



e (s 3 ) = 2 







slog(z )ds — s 3 / s \og(z )ds 



(i 



-C1S3 2 + C2S3 . 



(5.16) 



We see immediately that the constants must be zero to preserve analyticity near 
s 3 = 0. 

Following |10j we write both integrals in terms of z , determined as a function of 
s by equation (1.14 1: 



l = z 2 -72s 2 zg. 
One can then solve this equation for s as a function of zq : 

1 (z 2 -3) 



72z 3 



sds = 



144 zl 



dzn 



(5.17) 



(5.18) 



Thus we find: 



eo(S3)= "(^Fi) 



(-o 2 -l) 



z ° (z 2 -3) 

j — log(z)dz- 

z 4 



(^o 2 -i) 2 A 



z ° (^ 2 -l)(^-3) 

£ log(z)dz 



2;' 



■feMi° e (.,) + ^°- 1)2 f° +1) 

Zq d Zq 

^^^ 10g(Zo)+ 12^^ 



(*§-l) 2 

_1 1 (z -l)(z 2 -6z -3) 

-2 lo ew + i2 ^Tl) 

5.2. Example: g=l 

In the case of 3 = 1, we have that 



(5.19) 



(5.20) 



of manipulation using the first component of equation (4.6), we have 



where z±(s) is given by Proposition 4.2 as follows: from that proposition, after a bit 

(5.21) 



3 3 

fl(s,w) = W _1 Zi(s«7 1/2 ) =--8f h 0tW fo,ww - ^ s fo h O,w > 
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into which we subsitute t he s elf-similar scalings fo — wzq^sw 1 ^ 2 ) and h$ 
to find 



( see Theorem 



3.i 



expand the w-derivatives and then set w — 1 



-u - 



i s z o 



;sz 



-u 



-SUr 



(5.22) 



finally we use the algebraic relations (1.7| and (1.8), and the expression (5.18) to 
eliminate all but zq from the formula for z\ , and we have 



l (sg-l) 2 (zg+9)z 
11 j " 4 (z 2 -3)4 



Theorem |5.1| gives 
ei(«s) = 2 



i Hi(zo)ds- 



l Hi(zo)ds 



c lS 2 3 +c 2 . 



(5.23) 



(5.24) 



We make the change of variables, as before, to integrals with respect to zq using (5.18 ). 
We then have: 



(«g-l) f Z0 z 2 (z 



(z 2 



'^-Hx(z)dz- 



+C lS 2 + C 2 



(^ 2 -3) 
z(z 2 -l) 



H 1 (z)dz+ (5.25) 



A direct calculation of Hi(zq) using (5.17), (5.19) and (5.23) shows that the inte- 
grands of both integrals in (5.25) are regular at Zq = 1 (which corresponds to S3 = 0). 
Hence, the vanishing of ei(0), which follows from (1.17), implies that C 2 —0. To de- 
termine Ci we will need to appeal to the Taylor expansion of e\ and its combinatorial 
interpretation. 

Using the explicit form of H\ as a function of zq{s 3 ) one finds from (5.25) that 

(^o 2 -l) 



ei(s 3 ) : 



1 

24 
1 

24 



log 
log 



: lOg 

24 S 



72 

C x - 



1 

"48 
3 
2 



(5.26) 
(5.27) 
(5.28) 



where we have choosen C% — 3/2 so that the coefficient of s 2 in ei(ss) will be 3/2!, as 
can be checked by directly calculating the second derivative of (5.28) with respect to 
S3 using the differential relation in (5.18) and then evaluating at S3 = (equivalently 
zq = X). This value of the second order coefficient is required by the fact that there 
are three genus one maps with two vertices of valence 3. 

The procedure can be continued for as far as one wishes; the only real constraint 
is the ability to find explicit expressions for the z g needed. We state without the 
computations the formula derived for 62(53): 



1 (z 2 -l) 3 (4z 4 - 



-93 2 2 - 



261) 



960 



(^o 2 -3) 5 



(5.29) 



The formulas we have derived for eo, ex, and e 2 as functions of zq have much in common 
with those found in the case of even times in |10| . This suggests an extension of the 
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global result, proven for the case of even times in [7]); namely, we expect that for 
g>l, e g (s3) is a rational function of z 2 with singularities occurring only at z 2 = 3, 
where zq is related to S3 by 



1: 



■ 72s 3 Zq 



5.3. Taylor Coefficients of eo(s3) and ei(s 3 ) 

From formulas for e g in terms of zq it is a straightforward procedure to derive 
expressions for the Taylor coefficients using contour integrals. The trick as in [10] is 
to again change variables to integrals with respect to z . 

One represents the Taylor coefficients of eo(ss) as contour integrals and substitutes 
the expression (5.191 in terms of zq. The 2jth Taylor coefficient is given by: 

eo 



K. 



(0) 
2j 



(2i)! 



2iri 
1 1 
j 2iri 

1 1 



s 3 ~0 *3 

— s 

1 dz 3 



dS3 



d6 ° " 2 ^ 2 



(z 2 -3)(z 2 -2z-l 
6z(z + l) 2 

32i-l 2 3j-2 1 



72 



2: 



dz 



J 



2-ki 



z 3 i- 1 (z-l) 2 (z 2 -3)(z 2 -2z-l) 



dz 



(z 2 - iy+ 2 



(5.30) 
(5.31) 
(5.32) 
(5.33) 



where we have used an integration by parts to find line (5.31 ), and the notation 

I I J z~l 

indicates that the integral is over a small circle in the complex plane, oriented counter 
clockwise, and containing z — 1. We then note that 

(z - l) 2 (z 2 - 3) (z 2 - 2z - 1) = z 6 - 4z 5 + z 4 + 12z 3 - I3z + 3 . (5.34) 
We insert this into the integral together with the change of variables £ = z 2 : 
#2? _ 3 2 -''- 1 2 3 J- 2 




j r(|)r(3+i) 



(5.35) 



Likewise we can express the Taylor coefficients of the genus one expansion as 
contour integrals: 



K. 



32j-l 2 3j-2 

j2iri 



c (3 J+ i)/2 dc 

" (C-3) (C-i)^ 



(5.36) 
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